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There are two neutral B^-B^ meson systems, and B^^-B^ (generically denoted 

q = s,d), which exhibit particle-antiparticle mixing [1]. This mixing phenomenon is described 
in Ref. [2]. In the following, we adopt the notation introduced in Ref. [2] and assume CPT 
conservation throughout. In each system the light (L) and heavy (H) mass eigenstates, 

\B^,u)=p\B^,)±qm), (1) 

have a mass difference Am^ = mn—mi, > and a total decay width difference AF^ = Fh— Tl- In 
the absence of CP violation in the mixing, \q/p\ = 1, these differences are given by Arrig = 2|Mi2| 
and lAFgl = 2|ri2|, where M12 and F12 arc the off-diagonal elements of the mass and decay 
matrices [2]. The evolution of a pure \B^) or \B^) state at t = is given by 

|S^(t))=5+(t)|5^) + ^5-(i)l^°), (2) 

mit)) = g+{t)m) + ^^g.it)\B'^) , (3) 

which means that the flavor states remain unchanged (+) or oscillate into each other (— ) with 
time-dependent probabilities proportional to 



cosh ^ 2 ±cos(Amgt) 



(4) 



where Fo = (Fh + Fl)/2. In the absence of CP violation, the time-integrated mixing probability 
/ \g^{tf dt/{J \g-{t)f dt + J \g+{t)f dt) is given by 



xl + yl Amg ATq 



Standard Model predictions and phenomenology 

In the Standard Model, the transitions B^^B^ and B^—^Bg are due to the weak interaction. 
They are described, at the lowest order, by box diagrams involving two W bosons and two up- 
type quarks (see Fig. 1), as is the case for K^-K^ mixing. However, the long range interactions 
arising from intermediate virtual states are negligible for the neutral B meson systems, because 
the large B mass is off the region of hadronic resonances. The calculation of the dispersive and 
absorptive parts of the box diagrams yields the following predictions for the off-diagonal element 
of the mass and decay matrices [3], 
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Figure 1: Dominant box diagrams for the B^^B^ transitions (g = d or s). Similar diagrams 
exist where one or both t quarks are replaced with c or tt quarks. 



where Gj? is the Fermi constant, m^/- the W boson mass and mj the mass of quark i; m^^, /s^ 
and BBq are the B^ mass, weak decay constant and bag parameter, respectively. The known 
function 5*0 (xj) can be approximated very well by 0.784 x?-'^*^ [4] and Vij are the elements of the 
CKM matrix [5]. The QCD corrections •qs and -q'^ are of order unity. The only non negligible 
contributions to M12 are from box diagrams involving two top quarks. The phases of M12 and 
satisfy 

(8) 



mt 



implying that the mass eigenstates have mass and width differences of opposite signs. This 
means that, like in the K^-K^ system, the heavy state has a smaller decay width than that of 
the light state. Hence AF is expected to be negative in the Standard Model. 
Furthermore, the quantity 
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is small, and a power expansion of yields 
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Therefore, considering both Eqs. (8) and (9), the CP-violating parameter 
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is expected to be very smah: ~ 0(10"'^) for the 5^-5^ system and . 0(10"'') for the -B^-S^ 
system [6]. 

In the approximation of negligible CP violation in mixing, the ratio ATq/Aniq is equal to 
the small quantity IF12/M12I of Eq. (9); it is hence independent of CKM matrix elements, i.e. 
the same for the 52--B° and systems. It can be calculated with lattice QCD techniques; 

typical results are ~ 5x 10~^ with quoted uncertainties of ~ 30%. Given the current experimental 
knowledge (discussed below) on the mixing parameter Xg, 



Xd = 0.771 ± 0.012 
Xs > 20.6 at 95% CL 



(^2-^2 system 



B^ system) 



(12) 



the Standard Model thus predicts that AF^/r^ is very small (below 1%), but AF^/F^ consider- 
ably larger (~ 10%). These width differences are caused by the existence of final states to which 
both the S° and S° mesons can decay. Such decays involve b ccq quark-level transitions, 
which are Cabibbo-suppressed if q = d and Cabibbo-allowed if q 



s. 
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Experimental issues and methods for oscillation analyses 



Time-integrated measurements of B^-B mixing were published for the first time in 1987 by 
UAl [7] and ARGUS [8], and since then by many other experiments. These measurements 
are typically based on counting same-sign and opposite-sign lepton pairs from the semileptonic 
decay of the produced hh pairs. Such analyses cannot easily separate the contributions from the 
different 6-hadron species, therefore the clean environment of T(4S') machines (where only B^ 
and charged B^ mesons are produced) is in principle best suited to measure Xd- 

However, better sensitivity is obtained from time-dependent analyses aimed at the direct 
measurement of the oscillation frequencies Amcj and Am^ from the proper time distributions of 
B^j, or B^g candidates identified through their decay in (mostly) flavor-specific modes and suitably 
tagged as mixed or unmixed. (This is particularly true for the B^-B^ system where the large 
value of Xs implies maximal mixing, i.e. Xs — V^-) In such analyses the 5° or 5° mesons are 
either fully reconstructed, partially reconstructed from a charm meson, selected from a lepton 
with the characteristics of a b ^ i~ decay, or selected from a reconstructed displaced vertex. 
At high-energy colliders (LEP, SLC, Tevatron), the proper time t = is measured from the 
distance L between the production vertex and the B decay vertex, and from an estimate of the B 
momentum p. At asymmetric B factories (KEKB, PEP-II), producing e~^e~ —>■ T(4S) — > B^B^ 
events with a boost (= 0.425, 0.55), the proper time difference between the two B candidates 
is estimated as At ~ where Az is the spatial separation between the two B decay vertices 
along the boost direction. In all cases, the good resolution needed on the vertex positions is 
obtained with silicon detectors. 

The average statistical significance 5 of a or Bg oscillation signal can be approximated 
as [9] 

S « /sig (1 - 2ry) e-^^-'^*)'/^ , (13) 

where N and /gig are the number of candidates and the fraction of signal in the selected sample, ry 
is the total mistag probability, and at is the resolution on proper time (or proper time difference). 

The quantity S decreases very quickly as Am increases; this dependence is controlled by at, which 
is therefore a critical parameter for Amg analyses. At high-energy colliders, the proper time 
resolution at ~ "^^L © includes a constant contribution due to the decay length resolution 
a I, (typically 0.05-0.3 ps), and a term due to the relative momentum resolution a^jp (typically 
10-20% for partially reconstructed decays) which increases with proper time. At B factories, 
the boost of the B mesons is estimated from the known beam energies, and the term due to the 
spatial resolution dominates (typically 1-1.5 ps because of the much smaller B boost). 

In order to tag a B candidate as mixed or unmixed, it is necessary to determine its flavor 
both in the initial state and in the final state. The initial and final state mistag probabilities, r]i 
and r/j, degrade 5 by a total factor (1 — 2r/) = (1 — 2r/j)(l — 2r/j). In lepton-based analyses, the 
final state is tagged by the charge of the lepton from 6 — > £~ decays; the biggest contribution 
to 77/ is then due to 6 — ^ c — > decays. Alternatively, the charge of a reconstructed charm 
meson from or from i?^), or that of a kaon thought to come from a 6 ^ c ^ s 

decay [10], can be used. For fully inclusive analyses based on topological vertexing, final state 
tagging techniques include jet charge [11] and charge dipole [12] [13] methods. 

At high-energy colliders, the methods to tag the initial state (i.e. the state at production), 
can be divided in two groups: the ones that tag the initial charge of the h quark contained 
in the B candidate itself (same-side tag), and the ones that tag the initial charge of the other 
h quark produced in the event (opposite-side tag). On the same side, the charge of a track 
from the primary vertex is correlated with the production state of the B if that track is a 
decay product of a B** state or the first particle in the fragmentation chain [14] [15]. Jet- and 
vertex-charge techniques work on both sides and on the opposite side, respectively. Finally, the 
charge of a lepton from 6 — > ^~ or of a kaon from h c s can be used as opposite side 
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tags, keeping in mind that their performance is degraded due to integrated mixing. At SLC, the 
beam polarization produced a sizeable forward-backward asymmetry in the Z ^ bb decays and 
provided another very interesting and effective initial state tag based on the polar angle of the 
B candidate [12]. Initial state tags have also been combined to reach r/j ~ 26% at LEP [15] [16], 
or even 22% at SLD [12] with full efficiency. In the case ijf = 0, this corresponds to an effective 
tagging efficiency (defined as Q = e(l — 2r])^ where e is the tagging efficiency) in the range 
23 — 31%. The equivalent figure at CDF is ~ 3.5% for Tevatron Run I [17] (expected to reach 
^ 5% for Run II [18]), refiecting the fact that tagging is very challenging at hadron colliders. 

At B factories, the flavor of a B*^ meson at production cannot be determined, since the two 
neutral B mesons produced in a T{AS) decay evolve in a coherent P-wave state where they 
keep opposite flavors at any time. However, as soon as one of them decays, the other follows a 
time-evolution given by Eqs. (2) or (3), where t is replaced with At (which will take negative 
values half of the time). Hence, the "initial state" tag of a B can be taken as the final state 
tag of the other B. Effective tagging efficiencies Q = 28 — 29% are achieved by BABAR and 
Belle [19], using different techniques including h ^ f.^ and h ^ c ^ s tags. It is interesting to 
note that, in this case, mixing of this other B {i.e. the coherent mixing occurring before the 
first B decay) does not contribute to the mistag probability. 

In the absence of experimental evidence for a decay-width difference, oscillation analyses 
typically neglect AF in Eq. (4) and describe the data with the physics functions Fe~'"*(l it 
cos(Ami))/2 (high-energy colliders) or Fe~^l^*l(l ± cos(AmAi))/4 (asymmetric T(45) ma- 
chines). As can be seen from Eq. (4), a non-zero value of AF would effectively reduce the 
oscillation amplitude with a small time-dependent factor that would be very difficult to distin- 
guish from time resolution effects. Measurements of Arrid are usually extracted from the data 
using a maximum likelihood fit. No significant B^-B^ oscillations have been seen so far. To 
extract information useful to set lower limits on Arug, B^ analyses follow a method [9] in which 
a 5° oscillation amplitude A is measured as a function of a fixed test value of Amg, using a 
maximum likelihood fit based on the functions Fse~^''*(l it ^cos(Amst))/2. To a very good 
approximation, the statistical uncertainty on A is Gaussian and equal to 1/5 from Eq. (13). If 
AttIs = Atv^J^'^, one expects A = 1 within the total uncertainty o"^; however, if Am^ is (far) 
below its true value, a measurement consistent with ^ = is expected. A value of Amg can be 
excluded at 95% CL if ^ -|- 1.645 cr^ < 1. If Am^g^^ is very large, one expects ^ = 0, and all 
values of Am^ such that 1.645 (t^( Amg) < 1 are expected to be excluded at 95% CL. Because 
of the proper time resolution, the quantity (T^(Ams) is an increasing function of Amg and one 
therefore expects to be able to exclude individual Amg values up to Arrff^, where Am|^°®, 
called here the sensitivity of the analysis, is defined by 1.645 (T^(Am|^°^) = 1. 

mixing studies 

Many S^-B^ oscillations analyses have been published by the ALEPH [20], BABAR [21], 
Belle [22], CDF [14], DELPHI [13] [23], L3 [24] and OPAL [25] collaborations. Although a 
variety of different techniques have been used, the individual Anid results obtained at high- 
energy colliders have remarkably similar precision. Their average is compatible with the recent 
and more precise measurements from asymmetric B factories. The systematic uncertainties are 
not negligible; they are often dominated by sample composition, mistag probability, or 6-hadron 
lifetime contributions. Before being combined, the measurements are adjusted on the basis of 
a common set of input values, including the 6-hadron lifetimes and fractions published in this 
Review. Some measurements are statistically correlated. Systematic correlations arise both from 
common physics sources (fragmentation fractions, lifetimes, branching ratios of b hadrons), and 
from purely experimental or algorithmic effects (efficiency, resolution, tagging, background de- 
scription). Combining all published measurements [13,14,20-25] and accounting for all identified 
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correlations as described in Ref. [26] yields Am^ = 0.502 it 0.004(stat) ± 0.005(syst) ps ^. 

On the other hand, ARGUS and CLEO have published time-integrated measurements [27] 
[28] [29], which average to Xd = 0.182±0.015. Following Ref. [29], the width difference AF^ could 
in principle be extracted from the measured value of r,^ and the above averages for Am^ and Xd 
(see Eq. (5)), provided that AT^ has a negligible impact on the Am^ measurements. However, di- 
rect time-dependent studies yield stronger constraints: DELPHI published the result | AFdl/Fd < 
18% at 95% CL [13], while BABAR recently obtained -8.4% < sign(ReAcp)AFd/Fd < 6.8% at 
90% CL [30]. 

Assuming AF,^ = and no CP violation in mixing, and using the measured B'^ lifetime, the 
Anid and Xd results are combined to yield the world average 



Anid = 0.502 ± 0.007 ps~ 



(14) 



or, equivalently. 



Xd = 0.186 ±0.004. 



(15) 

Evidence for CP violation in mixing has been searched for, both with flavor-specific 
and inclusive 5^ decays, in samples where the initial flavor state is tagged. In the case of 
semileptonic (or other flavor-specific) decays, where the final state tag is also available, the 
following asymmetry [2] 



~ 1 - \q/p\ 



(16) 



has been measured, either in time-integrated analyses at CLEO [28] [29] [31] and CDF [32], or 
in time-dependent analyses at LEP [33] [34] [35] and BABAR [30] [36]. In the inclusive case, 
also investigated at LEP [34] [35] [37], no final state tag is used, and the asymmetry [38] 



Ar(^°(t)^all)-Ar(g(t)^all) 
N{B^,{t)^ all) + N{B^,{t)^ all) 



^ v4sL 



y sin(Amdt) 



sm 



/ Arudt 
V 2 



(17) 



must be measured as a function of the proper time to extract information on CP violation. In 

all cases asymmetries compatible with zero have been found, with a precision limited by the 
available statistics. A simple average of all published results for the B^ meson [29,31,33,35-37] 
yields .4sl = +0.002 it 0.013, or \q/p\d = 0.999 ± 0.006, a result which does not yet constrain 
the Standard Model. 

The AiUd result of Eq. (14) provides an estimate of 2|Mi2| and can be used, together with 
Eq. (6), to extract the magnitude of the CKM matrix clement Vtd within the Standard Model [39]. 
The main experimental uncertainties on the resulting estimate of \Vtd\ come from mt and Amd', 
however, the extraction is at present completely dominated by the uncertainty on the hadronic 
matrix element jBj_\/BB^ = 221 ± 28l!l22 MeV obtained from lattice QCD calculations [40]. 



Br^ mixing studies 

oscillations have been the subject of many studies from ALEPH [15] [41], CDF [42], 
DELPHI [13] [16] [43] [44], OPAL [45] and SLD [12] [46] [47]. No oscillation signal has been 
found so far. The most sensitive analyses appear to be the ones based on inclusive lepton 
samples at LEP. Because of their better proper time resolution, the small data samples analyzed 
inclusively at SLD, as well as the few fully reconstructed Bg decays at LEP, turn out to be also 
very useful to explore the high Am^ region. 

All results are limited by the available statistics. They can easily be combined, since all 
experiments provide measurements of the 5° oscillation amplitude. All published results [12] 
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[13] [16] [41] [42] [43] [45] [46] are averaged using the procedure of Ref. [26] to yield the combined 
amplitudes A shown in Fig. 2 as a function of Amg. The individual results have been adjusted 
to common physics inputs, and all known correlations have been accounted for; the sensitivities 
of the inclusive analyses, which depend directly through Eq. (13) on the assumed fraction fs of 
5° mesons in an unbiased sample of weakly-decaying b hadrons, have also been rescaled to a 
common average of fs = 0.107 ± 0.011. The combined sensitivity for 95% CL exclusion of Anis 
values is found to be 17.8 ps~^. All values of Amg below 14.4 ps~^ are excluded at 95% CL, 
which we express as 

Arus > UA ps"^ at 95% CL . 

The values between 14.4 and 21.8 ps~^ cannot be excluded, because the data is compatible with 
a signal in this region. However, no deviation from ^ = is seen in Fig. 2 that would indicate 
the observation of a signal. 

Some Arris analyses are still unpublished [44] [47] . Including these in the above combination 
would yield Amg > 14.6 ps~^ at 95% CL with a sensitivity of 19.3 ps~^. 

The information on IV^^I obtained, in the framework of the Standard Model, from the com- 
bined amplitude spectrum is hampered by the hadronic uncertainty, as in the case. However, 
several uncertainties cancel in the frequency ratio 

Ams ^ niB, ^2 
AiTid ruBa 

where ^ = (/b, \/BB, )/{ fBd ^/Bb^) = 1-15 ± 0.05lo;Jo is an SU(3) flavor-symmetry breaking 
factor obtained from lattice QCD calculations [40]. The CKM matrix can be constrained using 
the experimental results on Anid, Anig, \Vub/Vcb\, and sin(2/3) together with theoretical 
inputs and unitarity conditions [39] [48]. Given all measurements other than Ama and Arris, 
the constraint from our knowledge on the ratio Amg/ Arud is presently more effective in limiting 
the position of the apex of the CKM unitarity triangle than the one obtained from the Am^ 
measurements alone, due to the reduced hadronic uncertainty in Eq. (18). We note also that it 
would be difficult for the Standard Model to accommodate values of Arris above ~ 25 ps^^ [48]. 

Information on AFs can be obtained by studying the proper time distribution of untagged 
data samples enriched in 5° mesons [49]. In the case of an inclusive B° selection [50] or a 
semileptonic B^^ decay selection [16] [51], both the short- and long-lived components are present, 
and the proper time distribution is a superposition of two exponentials with decay constants 
Fs ± AFs/2. In principle, this provides sensitivity to both Vg and (AFs/Fs)^. Ignoring AFs 
and fitting for a single exponential leads to an estimate of Fs with a relative bias proportional 
to (AFs/Fs)^. An alternative approach, which is directly sensitive to first order in AFg/Fs, 
is to determine the lifetime of candidates decaying to CP cigcnstates; measurements exist 
for Bl J/^tj) [52] and B^ [53], which are mostly CP-even states [54]. An 

estimate of AFs/Fs has also been obtained directly from a measurement of the 5° D^f''^ D^s^~ 
branching ratio [53], under the assumption that these decays account for all the CP-even final 
states (however, no systematic uncertainty due to this assumption is given, so the average quoted 
below will not include this estimate). 

Present data is not precise enough to efficiently constrain both Fs and AFs/Fs; since the 
Ps and P° lifetimes are predicted to be equal within less than a percent [55], an expectation 
compatible with the current experimental data [56], the constraint Fs = F^ can also be used 
to improve the extraction of AFs/Fs. Applying the combination procedure of Ref. [26] on the 
published results [16] [51-53, 57] yields 

|AFs|/Fs < 0.54 at 95% CL (19) 



Vtd 



(18) 
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Figure 2: Combined measurements of the oscillation amplitude as a function of Auig, in- 
cluding all results published by November 2003. The measurements are dominated by statistical 
uncertainties. Neighboring points are statistically correlated. 



without external constraint, or 

|Ar^|/r, < 0.29 at 95% CL (20) 

when constraining l/Tg to the measured lifetime. These results are not yet precise enough 
to test Standard Model predictions. 

Average 6-hadron mixing and 6-hadron production fractions at 
high energy 

Let fu, fd: fs and /baryon be the fractions of Bu, B^, B^ and 6-baryon composing an unbiased 
sample of weakly-decaying b hadrons produced in high-energy colliders. LEP experiments have 
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measured fs x BR{B^ Dj(.+ViX) [58], BR(6 ^ A^) x BR(A0 ^ K+t-V^X) [59] and BR(6 ^ 
S^) X BR(E;^ E~£~i'(X) [60] from partially reconstructed final states including a lepton, 
/baryon from protons identified in b events [61], and the production rate of charged b hadrons 
[62]. The various fe-hadron fractions have also been measured at CDF from electron-charm 
final states [63]. All these published results have been combined following the procedure and 
assumptions described in Ref. [26], to yield /„ = = (40.3 ± 1.1)%, = (8.8 ± 2.1)% and 
/baryon = (10.7 it 1.8)% uudcr the constraints 

fu = fd and fu + fd + fs + /baryon = 1 ■ (21) 

Time-integrated mixing analyses performed with lepton pairs from bb events produced at 
high-energy colliders measure the quantity 

X = fdXd + f'sXs, (22) 

where /^ and f'g are the fractions of and hadrons in a sample of semileptonic 6-hadron 
decays. Assuming that all b hadrons have the same semileptonic decay width implies = 
fq/{rqTh) {q = s,d), where is the average 6-hadron lifetime. Hence x measurements can be 
used to improve our knowledge on the fractions /„, fd, fs and /baryon- 

Combining the above estimates of these fractions with the average x = 0.1257 ± 0.0042 
(published in this Review), Xd from Eq. (15) and Xs = 1/2 yields, under the constraints of 
Eq. (21), 

fu = fd = (39.7 ± 1.0)% , (23) 
/, = (10.7 ±1.1)%, (24) 

/baryon = (9.9 ±1.7)%, (25) 

showing that mixing information substantially reduces the uncertainty on fg. These results and 
the averages quoted in Eqs. (14) and (15) for Xd and Arrid have been obtained in a consistent 
way by the B oscillations working group [26], taking into account the fact that many individual 
measurements of Am^ depend on the assumed values for the 6-hadron fractions. 



Summary and prospects 

B^~B^ mixing has been and still is a field of intense study. The mass difference in the B^-B^ 
system is very well measured (with an accuracy of 1.3%) but, despite an impressive theoretical 
effort, the hadronic uncertainty still limits the precision of the extracted estimate of \Vtd\- The 
mass difference in the B^-B^ system is much larger and still unmeasured. However, the current 
experimental lower limit on Am^ already provides, together with Anid, a significant constraint 
on the CKM matrix within the Standard Model. No strong experimental evidence exists yet for 
the rather large decay width difference expected in the B^—B^ system. It is interesting to recall 
that the ratio ATg/ Arris does not depend on CKM matrix elements in the Standard Model (see 
Eq. (9)), and that a measurement of either Arris or AT^ could be turned into a Standard Model 
prediction of the other one. 

In the near future, the most promising prospects for B^ mixing are from Run II at the 
Tevatron, where both Arris and AFg are expected to be measured with fully reconstructed B^ 
decays. The CDF and DO collaborations expect to be able to observe B^ oscillations with 
2 — 3 fb""*^ of data, if Arris is consistent with the current Standard Model prediction [18]. Should 
this not be the case, then the discovery of B^ oscillations will most likely be made at CERN's 
Large Hadron Collider scheduled to come into operation in 2007, where the LHCb collaboration 
claims to have the potential to cover a Arris range up to ~ 68 ps-i after 2 fb"^ of data (lO'^ s) 
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have been analyzed [64] . The BTeV experiment at Fermilab should have a comparable sensitivity 
to Ams and is expected to turn on in 2009 [65] . 

CP violation in B mixing, which has not been seen yet, as well as the phases involved in B 
mixing, will be further investigated with the large statistics that will become available at the B 
factories, the Tevatron and the LHC. 

B mixing may not have delivered all its secrets yet, because it is one of the phenomena 
where new physics might very well reveal itself (for example new particles involved in the box 
diagrams). Theoretical calculations in lattice QCD are becoming more reliable and further 
progress in reducing hadronic uncertainties is expected. In the long term, a stringent check 
of the consistency, within the Standard Model, of the S° and S° mixing measurements with 
all other measured observables in B physics (including CP asymmetries in B decays) will be 
possible, allowing to place limits on new physics or, better, discover new physics. 
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